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Abstract. In an XMM-Newton raster observation of the bright Local Group spiral galaxy M 33 we study the population of 
X-ray sources (X-ray binaries, supernova remnants) down to a 0.2-4.5 keV luminosity of 10^' erg s"' - more than a factor 
of 10 deeper than earlier ROSAT observations. EPIC hardness ratios and optical and radio information are used to distinguish 
between different source classes. The survey detects 408 sources in an area of 0.80 square degree. We correlate these newly 
detected sources with earlier M 33 X-ray catalogues and information from optical, infra-red and radio wavelengths. As M 33 
sources we detect 21 supernova remnants (SNR) and 23 SNR candidates, 5 super-soft sources, and 2 X-ray binaries (XRBs). 
There are 267 sources classified as hard, which may either be XRBs or Crab-like SNRs in M 33 or background AGN. The 44 
confirmed and candidate SNRs more than double the number of X-ray detected SNRs in M 33. 16 of these are proposed as 
SNR candidates from the X-ray data for the first time. On the other hand, there are several sources not connected to M 33: 
five foreground stars, 30 foreground star candidates, 12 active galactic nucleus candidates, one background galaxy and one 
background galaxy candidate. Extrapolating from deep field observations we would expect 175 to 210 background sources in 
this field. This indicates that about half of the sources detected are sources within M 33. 
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1. Introduction 

The Local Group Sc galaxy M 33 is located at a distance of 795 
kpc (van den Bergh 1991, i.e. 1" corresponds to 3.9 pc and the 
flux to luminosity conversion factor is 7.6x10'*'^ cm^), is seen 
under a relatively low inclination of 56° (Zaritsky et al. 1989) 
and the optical extent can be approximated by an inclination- 
corrected D25 ellipse with large diameter of 64 .'4 and axes ra- 
tio of 1.66 (de Vaucouleurs et al. 1991; Tully 1988). With its 
moderate Galactic foreground absorption (A^h= 6 x 10^" cm"^. 
Stark et al. 1992) M 33 is ideally suited to study the X- 
ray source population and diffuse emission in a nearby spiral 
galaxy. The Einstein X-ray Observatory detected diffuse emis- 
sion from hot gas in M 33 and 17 unresolved sources (Long 
et al. 1981; Markert & Rallis 1983; Trinchieri et al. 1988). First 
ROSAT HRI and PSPC observations revealed 57 sources and 
confirmed the detection of diffuse X-ray emission which may 
trace the spiral arms within a 10' radius around the nucleus 
(Schulman & Bregman 1995; Long et al. 1996). Combining all 
archival ROSAT observations of the field, Haberl & Pietsch 
(2001, hereafter HPOl) found 184 X-ray sources within a 50' 
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radius around the nucleus, identified some of the sources by 
correlations with previous X-ray, optical and radio catalogues, 
and in addition classified sources according to their X-ray 
properties. They found candidates for super-soft X-ray sources 
(SSS), X-ray binaries (XRBs), supernova remnants (SNRs), 
foreground stars and active galactic nuclei (AGN). 

To follow up on these findings, we proposed as part of the 
telescope scientist guaranteed time (GT) program a deep ho- 
mogeneous XMM-Newton survey of M 33 with a sensitivity of 
lO"*^ erg s"' in the 0.5-10 keV band, a factor of ten deeper than 
previous surveys. The survey consisted of a raster of 15 point- 
ings of about 10 ks each. The directions were selected in a way 
that each area within the M 33 optical D25 extent was covered 
by the XMM-Newton EPIC detectors at least three times with 
the medium filter. The survey was complemented by an addi- 
tional A02 program, initiated to fill observation gaps due to 
high detector background, or, because the thick filter was used 
during the GT observations. First results of this survey were 
presented in Pietsch et al. (2003a). 

Two M 33 sources are known for their outstanding X- 
ray properties (Peres et al. 1989). The brightest source (X- 
8 in the nomenclature of Long et al. 1981, luminosity of 



10 



,39-40 



erg s ) is the most luminous X-ray source in the Local 



Group of galaxies and coincides with the optical center of 
M 33. Its time variability (Dubus et al. 1997), its point-hke 
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Table 1. XMM-Newton M 33 observation log (proposal numbers 010264 and 014198). 



Field Obs.id. Obs. dates Pointing direction EPICPN EPICMOSl EPICM0S2 



V133 








RA/DEC (J2000) 


Mode+/Filter 




Mode+/Filter 




Mode+/Filter 












(4) 


(5) 




\i > 


\°) 








1 


0102640101 


2000-08-04 


1 


33:51.0 


30:39:37 


FF/medium 


6.29 


PW2/medium 


6.23 


PW2/medium 


6.22 


1 


0141980501 


2003-01-22 


1 


33:51.0 


30:39:37 


FF/medium 


1.23 


FF/medium 


1.51 


FF/medium 


1.51 


1 


0141980801 


2003-02-12 


1 


33:51.0 


30:39:37 


FF/medium 


7.25 


FF/medium 


8.86 


FF/medium 


8.87 


2 


0102640201 


2000-08-04 


1 


34:40.0 


30:57:48 


EFF/medium 


11.40 


FF/medium 


15.24 


FF/medium 


15.24 


3 


0102640301 


2000-08-07 


1 


33:32.0 


30:52:13 


EFF/medium 


3.80 


FF/medium 


4.58 


FF/thin 


4.58 


3 


0141980401 


2003-01-24 


1 


33:32.0 


30:52:13 


FF/medium 


0.00 


FF/medium 


0.00 


FF/medium 


0.00 


4 


0102640401 


2000-08-02 


1 


32:51.0 


30:36:49 


EFF/thick 


8.66 


FF/thick 


12.99 


FF/thick 


12.99 


4 


0141980601 


2003-01-23 


1 


32:51.0 


30:36:49 


FF/medium 


11.18 


FF/medium 


12.56 


FF/medium 


12.56 


5 


0102640501 


2001-07-05 


1 


33:02.0 


30:21:24 


FF/medium 


9.26 


FF/medium 


11.76 


FF/medium 


11.76 


6 


0102640601 


2001-07-05 


1 


34:08.0 


30:46:06 


FF/medium 


2.86 


FF/medium 


3.42 


FF/medium 


3.42 


6 


0141980301 


2003-07-25 


1 


34:08.0 


30:46:06 


FF/medium 


4.12 


FF/medium 


4.82 


FF/medium 


4.82 


7 


0102640701 


2001-07-05 


1 


34:10.0 


30:27:00 


FF/medium 


6.97 


FF/medium 


10.96 


FF/medium 


10.95 


8 


0102640801 


2001-07-07 


1 


34:51.0 


30:42:22 


FF/medium 


- 


FF/medium 


1.19 


FF/medium 


1.19 


8 


0102642001 


2001-08-15 


1 


34:51.0 


30:42:22 


FF/medium 


8.65 


FF/medium 


11.15 


FF/medium 


11.15 


9 


0102640901 


2001-07-07 


1 


34:04.0 


30:57:25 


FF/medium 


2.42 


FF/medium 


3.15 


FF/medium 


3.15 


9 


0141980201 


2003-07-11 


1 


34:04.0 


30:57:25 


FF/medium 


3.74 


FF/medium 


4.02 


FF/medium 


4.02 


10 


0102641001 


2001-07-08 


1 


33:07.0 


30:45:02 


IT/medium 


1.25 


FF/medium 


8.33 


FF/medium 


8.37 


10 


0141980101 


2003-07-11 


1 


33:07.0 


30:45:02 


FF/medium 


5.45 


FF/medium 


5.75 


FF/medium 


5.77 


11 


0102641101 


2001-07-08 


1 


32:46.0 


30:28:19 


FF/medium 


8.10 


FF/medium 


10.36 


FF/medium 


10.36 


12 


0102641201 


2000-08-02 


1 


33:38.0 


30:21:49 


EIT/thick 


11.34 


FF/thick 


3.67 


FF/thick 


3.66 


12 


0141980701 


2003-01-24 


1 


33:38.0 


30:21:49 


IT/medium 


4.40 


FF/medium 


5.72 


FF/medium 


5.73 


13 


0102642101 


2002-01-25 


1 


34:34.0 


30:34:11 


IT/medium 


10.00 


FF/medium 


12.27 


FF/medium 


12.27 


14 


0102642201 


2002-01-25 


1 


34:56.0 


30:50:52 


FF/medium 


11.60 


FF/mcdium 


13.87 


FF/mcdium 


13.87 


15 


0102642301 


2002-01-27 


1 


33:33.0 


30:33:07 


FF/medium 


10.00 


FF/medium 


12.27 


FF/medium 


12.27 



Notes: 

* : FF: full frame, EFF: extended full frame. PW2: partial window 2 

t : Exposure time in units of ks after screening for high background (see text) 



nature seen by ROSAT HRI (Pietsch & Haberl 2000) and 
Chandra (Dubus & Rutledge 2002) and its X-ray spectrum best 
described by an absorbed power law plus disk blackbody model 
(e.g. Ehle et al. 2001; La Parola et al. 2003; Foschini et al. 
2003), point towards a black hole XRB. A possible periodic- 
ity of 106 days was not confirmed in later observations (see 
Parmar et al. 2001). The second source (X-7) is an ecUpsing 
high mass XRB (HMXB) with a binary period of 3.45 d and 
possible 0.31 s pulsations, discovered in ROSAT PSPC and 
HRI observations (Schubnan et al. 1993; Dubus et al. 1997, 
1999; Larson & Schulman 1997). Pietsch et al. (2004) ana- 
lyzed XMM-Newton and Chandra data and determined im- 
proved binary parameters, however, they could not confirm the 
proposed pulsations. In a special analysis of DIRECT data, they 
identified an 071 star of 18.89 mag in V, as the optical coun- 
terpart, which shows the ellipsoidal heating hght curve of a 
HMXB with the X-7 binary ephemeris. X-7 was the most dis- 
tant eclipsing XRB until the detection of another such source 
in XMM-Newton and Chandra observations of the starburst 
galaxy NGC 253, which is located at more than three times the 
distance of M 33 (Pietsch et al. 2003b). 

Here we present merged medium and thin filter images for 
the three EPIC instruments, in five energy bands, using only 



Fig. 1. XMM-Newton EPIC "PN equivalent" exposure map combin- 
ing all low background observation times of the fields of the M 33 
raster. The outer contour indicates ks, gray-scale steps every 10 ks. 
The maximum PN equivalent exposure (see text) is 89 ks in the cen- 
ter area of M 33. The optical extent of M 33 is marked by the white 
inclination-corrected D25 ellipse. The center of the fields of our raster 
(see Table 1) are indicated and numbered. 

times of low background from our XMM-Newton survey of 
M 33, as well as a source catalogue and source population study 
based on these images. An analysis of the individual pointings 
with less stringent background rejection and using all filters 
is in progress, and will result in a source catalogue with im- 
proved positions (taking care of individual pointing offsets), 
that will also provide information on source variability and 
spectral characteristics. 

2. Observations and data analysis 

Table 1 summarizes the XMM-Newton (Jansen et al. 2001) 
EPIC (Struder et al. 2001; Turner et al. 2001) observations. For 
each observation we give the field number of our M 33 raster 
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Table 2. Count rate to energy conversion factor for different spectral 
models for EPIC instruments in energy band 1 to 5 (B1-B5). Models 
are a power law with photon index of 1.7 (PL), a thin thermal model 
with temperature of 1 keV (TT), and a black body with a temperature 
of 30 eV (BB). All models assume galactic foreground absorption of 
6x10^° cm-2 (Stark et al. 1992). 



Detector 


Spec 


Bl 


B2 


B3 


B4 


B5 








(10-1^ 


erg cm"^ 


ct-i) 




EPIC PN 


PL 


1.03 


1.10 


1.71 


4.9 


21. 




TT 


1.03 


1.14 


1.55 


4.0 






BB 


1.09 


0.97 


1.30 






EPIC MOS 


PL 


6.2 


4.9 


4.8 


12.9 


66. 




TT 


6.0 


4.5 


4.6 


11.0 






BB 


10.1 


4.8 


4.2 







(col. 1), the observation identification (2), date (3), pointing di- 
rection (4,5), as well as instrument mode, filter and exposure 
time after screening for high background for EPIC PN (6,7), 
MOSl (8,9), and M0S2 (10,1 1). Observation 0102640801 was 
affected by a satellite "slew failure" and only the EPIC MOS 
instruments got some exposure time. 

In the XMM-Newton observations the EPIC PN and MOS 
instruments were mostly operated in the fuU frame mode re- 
sulting in a time resolution of 73.4 ms and 2.6 s, respectively. 
Only for four early observations the PN detector was oper- 
ated in the extended full frame mode (time resolution 200 ms). 
During observation 0102640101 the MOS detectors were op- 
erated in the small window mode (0.3 s time resolution for the 
iimer CCDs to avoid pile-up for the bright source X-8). The 
medium filter was in front of the EPIC cameras in all but the 
first two observations which were performed with the thick fil- 
ter (0102640401, 0102641201) and observation 0102640301 
for which the MOS2 detector was operated with the thin filter. 
For creating the merged images and for source detection, we 
only used medium and thin filter observations. 

The data analysis was performed using tools in the XMM- 
Newton Science Analysis System (SAS) v5.4.1 and some later 
versions from the development area as specially mentioned, 
EXSAS/MIDAS 1.2/1.4, and FTOOLS v5.2 software pack- 
ages, the imaging application DS9 v2.3bl and the spectral anal- 
ysis software XSPEC vl 1.2. 

For MOS, we used SAS task emchain with task emevertt 
v7.10 that removes flickering pixels. We carefully screened 
the event files for bad CCD pixels, remaining after the stan- 
dard processing. To create a homogeneous combined image 
with similar background level for all fields we had to care- 
fully screen the data for high background times. In most ob- 
servations the background light curves at all energies follow 
the high energy (7-15 keV) background light curves provided 
by the SAS tasks epchain and emchain and allowed us to effi- 
ciently screen the data. In some observations, however, the low 
energy background (below 1 keV) was higher during high en- 
ergy selected good time intervals (GTIs) and had to be used for 
additional screening. GTIs were determined from the higher 



signal-to-noise PN light curves and also used for the MOS 
cameras. Outside the PN time coverage, GTIs were determined 
from the combined MOS light curves. The corresponding low 
background times for the individual observations are listed for 
the PN and MOS cameras in Table 1 . The combined "PN equiv- 
alent" vignetted EPIC exposure (1-2 keV) is significantly above 
20 ks within most of the optical extent of M 33 (indicated by 
the inclination-corrected D25 ellipse in Fig. 1). Only in the NW 
and S the exposure does not reach this limit. For the "PN equiv- 
alent" EPIC exposure map we added MOSl and MOS2 expo- 
sure maps (scaled by a factor of 0.36 as expected for a 5 keV 
thermal bremsstrahlung spectrum with an absorption of 6x10^" 
cm"^, as expected for XRBs in M 33) to the PN exposure. 

For the analysis we used five energy bands: (0.2-0.5) keV, 
(0.5-1.0) keV, (1.0-2.0) keV, (2.0-4.5) keV, and (4.5-12) keV 
as band 1 to 5. We intentionally split the 0.5-2.0 keV band used 
in the IXMM XMM-Newton Serendipitous Source Catalogue^ 
to get, on average, a more homogeneous distribution of the 
source counts to the energy bands, which leads to a better 
spread of the hardness ratio values and allows effective source 
classification (see Sect. 6). 

For PN we selected only "singles" (PATTERN=0) in band 
1, for the other bands "singles and doubles" (PATTERN<4). 
For MOS we used "singles" to "quadruples" (PATTERN<12). 
To avoid background variabiUty over the PN images we omit- 
ted the energy range (7.2-9.2) keV in band 5 where strong flu- 
orescence lines cause higher background in the outer detector 
area (Freyberg et al. 2004). To convert source count rates in the 
individual bands to fluxes we calculated count rate to energy 
conversion factors (ECF) for PN and MOS medium filter ob- 
servations using W3PIMMS and assuming the same spectrum 
as for the first XMM-Newton source catalogue, i.e. a power law 
spectrum with photon index 1.7 absorbed by the galactic fore- 
ground column of 6x10^^* cm"^ (Table 2). The assumed spec- 
trum represents hard sources like XRBs or AGN. We also give 
ECF values for a typical SNR (absorbed 1 keV thin thermal) 
and SSS (absorbed 30 eV black body) in M 33. Table 2 demon- 
strates that for EPIC PN, ECF values only vary by about 20% 
for the different spectra in bands 1 to 4. The same is true for 
EPIC MOS, with the exception of the black body ECF in the 
0.2-0.5 keV band where the MOS sensitivity is much lower. 
For both PN and MOS, the sensitivity is strongly reduced at 
energies above 4.5 keV. On the other hand, background count 
rates are similar in the hard band. Therefore, only bright hard 
sources are detected in this band. 

To classify the source spectra we computed four hardness 
ratios from the source count rates. These hardness ratios are de- 
fined as HRi - (B,+i - B,) (B,+ i + B,) , for i = 1 to 4, where B, 
denotes count rates in band i, as defined above. In the standard 
source detection products, hardness ratios and count rates in in- 
dividual energy bands are not combined for all instruments. To 
improve the statistics, we have combined the count rates (and 
fluxes) from all three EPIC instruments and have then derived 
hardness ratios for them. We are aware that these products have 
some jitter in their meaning, as the relative integration times for 

' Prepared by the XMM-Newton Survey Science Centre 
Consortium (http : //xitimssc-www. star . le . ac. uk/) 
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individual sources in the EPIC instruments differ slightly. Due 
to the better signal in the combined products, they are still very 
valuable for source classification. For investigations needing 
accurate calibration, values derived for each instrument sepa- 
rately should be used. 

As discussed above, band 5 suffers from the combined ef- 
fect of lower sensitivity in both PN and MOS, and a large con- 
tamination from the hard background, reducing the significance 
of any detection due to the much higher noise. We therefore 
give total count rates and errors derived from all instruments 
in the "XID" band that comprises bands 1^, and derive fluxes 
in the corresponding 0.2-4.5 keV range. A larger band would 
artificially increase the overall noise of all but extremely hard 
sources. For very soft sources, or for sources detected only in 
one of the instruments (e.g. soft sources in PN), even with this 
choice the quoted errors are overestimated, however the num- 
ber of sources of this nature is small. 

For PN, MOSl and M0S2 we created images, background 
images, exposure maps (without and with vignetting correc- 
tion) and masked them for the acceptable detector area in each 
of the five energy bands. For PN, the background maps contain 
the contribution from the "out of time (OOT)" events (param- 
eter withootset=true in task esplinemap). The images of 
the individual observations were then merged for source detec- 
tion into images covering the full M 33 field. The maximum 
vignetted exposure was 44.3 ks, 58.6 ks and 58.7 ks in the 1- 
2 keV band for the PN, MOSl and MOS2 cameras, respec- 
tively, and was nearly the same for bands 1 to 4. In the hard- 
est band, the maximum exposure is smaller by 30% due to the 
larger vignetting efi'ects. To allow us an easy merging of im- 
ages, we calculated for the events of all observations projected 
sky coordinates (X,Y) with respect to the center reference po- 
sition RA=0l''33™5Oi4, DEC=+30°39'36"(J2000). At the bor- 
ders, due to the diff'erent field of view of the EPIC instruments, 
the merged images for PN, MOS 1 and MOS2 do not fully over- 
lap. Nevertheless, we used the fuU area for source detection. 

3. Source catalogue 

We searched for sources in these merged images using simul- 
taneously 5x3 images (5 energy bands and PN, MOSl and 
M0S2 camera). A preliminary source list created with the task 
boxdetect with a low likelihood threshold was used as start- 
ing point for the task emldetect v4 . 3 1. As photons collected 
with different off-axis angles contribute to the source counts in 
the merged image, we prepared a specific point spread function 
(PSF) calibration file. It assumed a constant PSF over the total 
field of view that was degraded as if 6' off'-axis. This is the most 
degraded but still "circular" PSF in the EPIC calibration files 
and a first approximation for the PSF when source counts are 
collected at many different off-axis angles as in the M 33 raster 
observation. To resolve sources that overlap due to the PSF, we 
used parameter multisourcef it=4. 

A detection probabiUty above 4cr is equivalent to a single 
band limiting detection UkeUhood of 10. This threshold has 
been successfully used by HPOl for the ROSAT PSPC detec- 
tions (see Sect. 5). In the case of simultaneous detection with 
a large number of images this equivalence is no longer valid 



due to the small number of counts in the individual images 
contributing to a source at the detection threshold. We have 
established for the task emldetect that for the case of simul- 
taneously using 15 images, a detection probability above 4cr 
corresponds to a combined UkeUhood of 7. This threshold was 
used to select sources for the catalogue. The detection proce- 
dure produced acceptable results as can be checked from the 
overlay of the detected sources onto the smoothed X-ray im- 
ages (see Sect. 4). 

Two emission peaks were resolved into several sources with 
no sign of multipUcity from the smoothed images and are better 
described as extended. In a detection run including extent de- 
termination we found that these two sources (#125 and #299) 
can be characterized by a Gaussian with FWHM of 18'.'0 and 
16'.'2, respectively, and the extent UkeUhood is high. In total we 
detected 408 sources in the field. 

The source parameters are summarized in Tables 3 (EPIC 
combined products) and 4 (products for EPIC PN, MOS 1 and 
MOS2, separately). For both tables we only show the first five 
and the last row in the paper. The fuU tables are available in 
electronic form at the CDS (Editors should add URL address 
here...). 

Table 3 gives the source number (col. 1), position (2 and 

3) with Icr uncertainty radius (4), likelihood of existence (5), 
integrated PN, MOSl and MOS2 count rate (6) and flux (7) in 
the 0.2^.5 keV band, and hardness ratios (8-1 1). Hardness ra- 
tios are calculated only for sources for which at least one of the 
two band count rates has a significance greater than 2cr. Errors 
are the properly combined statistical errors in each band and 
can extend beyond the range of allowed values of hardness ra- 
tios as defined previously (-1.0 to 1.0). The EPIC instruments 
contributing to the source detection, are indicated in the "Val" 
parameter (col. 12, first character for PN, second MOSl, third 
MOS2) as "T", if inside the field of view (FOV), or "F", if out- 
side of FOV. There are only nine sources at the periphery of 
the FOV where only part of the EPIC instruments contribute. 
The positional error does not include intrinsic systematic er- 
rors which, according to the Users Guide of the IXMM XMM- 
Newton Serendipitous Source Catalogue, amount to r.'5 and 
should be quadratically added to the statistical errors in col. 
4. In the remaining columns of Table 3, we give cross corre- 
lation information, which is further described in Sect. 5. We 
only want to mention here, that we used the foreground stars 
and candidates, to verify the assumed source position errors. 
All 35 foreground stars and candidates are located within the 
3cr statistical plus systematic positional error given above. 

The faintest sources detected have a flux of 1.0x10"'^ 
ergcm^^s \ the brightest source (X-8) 1.1x10"'^ 
erg cm"^ s"' in the 0.2^.5 keV band. This corresponds 
to an absorbed luminosity range in M 33 of 7.6x10^^ erg s~' 
to 8.4x10^^ ergs"'. 

Table 4 repeats source number and Val parameter of Table 3 
(cols. 1 and 2). It then gives for EPIC PN exposure (3), source 
existence likelihood (4), count rate (5) and flux (6) in the 0.2- 
4.5 keV band, and hardness ratios (7-10). Columns 11 to 18 
and 19 to 26 give the same information corresponding to cols. 
3 to 10, but now for the EPIC MOSl and MOS2 instruments. 
Hardness ratios were again screened as in Table 3. From the 
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Table 3. X-ray source catalogue of the XMM-Newton EPIC M 33 raster observation. 



5 



Src 


RA(J2000) 


DEC{J2000) 






Count rate* 


Flux' 


HRl 


HR2 


HR3 


HR4 


Val 




(hmmss.ss) 


{+dd mm ss.s) 


(") 




(cts-') 


(erg cm"- s"') 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


1 31 42.93 


-1-30 23 29.5 


1.78 


3.41e+01 


1.18e-02±2.2e-03 


1.55e-14±4.0e-15 


0.44±0.21 


-0.39+0.21 


-0.84±0.34 


*± • 


TFF 


2 


1 31 46.89 


-1-30 22 01.3 


1.03 


1.43e+02 


3.86e-02±4.1e-03 


4.17e-14±5.7e-15 


0.44±0.11 


-0.17±0.12 


-0.75±0.14 


•± • 


TTF 


3 


1 31 54.90 


+30 29 54.5 


1.86 


2.48e+01 


7.15e-03±1.6e-03 


6.33e-15±1.7e-15 


0.37±0.21 


-0.64+0.26 




* 


TTT 


4 


1 31 55.12 


+30 28 56.9 


1.62 


1.49e+01 


6.37e-03±1.7e-03 


7.72e-15±2.7e-15 




0.31+0.31 


-0.11+0.29 


0.20+0.35 


TTT 


5 


1 31 59.18 


+30 35 06.9 


1.60 


1.58e+01 


5.00e-03±1.0e-03 


6.96e-15±1.7e-15 


0.38+0.27 


-0.37+0.27 


0.37+0.28 


-1.00+0.90 


TTT 


408 


1 36 01.07 


+30 51 28.0 


1.10 


l.lOe+02 


2.52e-02±2.7e-03 


2.72e-14±4.1e-15 


0.82+0.13 


0.06+0.11 


-0.32+0.14 


-0.20+0.31 


TTT 



XID USNO-Bl Mul A B2 R2 12 108(177) '^'"^^ Remarks 



name (") (mag) (mag) (mag) 



(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 




1203-0020312 


1 


3.1 


15.8 


13.5 


11.3 


-2.6 


fgStar 


WUMaVar(SIM) 


LCB2(Var(flare)),HP8 


2293-01403-1 


1 


2.7 


11.1 


10.8 


10.6 


-3.6 


fgStar 


betaLyraeEB(SIM) 




2293-01131-1 


1 


2.7 


12.4 


11.3 


10.9 


-4.1 


< fgStar > 


<G8> 






* 


* 


* 


* 


* 


* 


< hard > 










* 


* 


* 




* 












* 


* 


* 


* 


* 


< hard > 





Notes: 



in XID band (0.2^.5 keV) 
References: 

in XID: X-n: Long et al. (1981), Trinchieri et al. (1988), SBn: Schuhnan & Bregman (1995), LCBn: Long et al. (1996), DCL97: Dubus et al. (1997), DCL99: Dubus et al. (1999), HPn: 
Haberl & Pietsch (2001) 

in remarks: SM: SIMBAD database, NED: NASA Extragalactic Database, NVSS: Condon et al. (1998), CS82: Christian & Schommer (1982), [MC83]25: Massey & Conti (1983), 
IFM*: Ivanov et al. (1993), GKLn: Gordon et al. (1998), GDKn: Gordon et al. (1999), IV74: Israel & van der Kruit (1974), VBW78: von Kap-Herr et al. (1978), M98: Massey (1998), 
MCM2001: Magrini et al. (2001), MBH96: Massey et al. (1996), PMM2004: Pietsch et al. (2004), GM94: Giovagnoh & Mould (1994), FSZ: Fabrika et al. (1997), D33: Macri et al. 
(2001),Mochejska et al. (2001), RW93: Regan & Wilson (1993), VHK: Hubble (1926), CGC97: Crampton et al. (1997), BCLMP: Boixlesteix et al. (1974) 



comparison of the hardness ratios in Table 3 and 4 it is clear that 
combining the instrument count rate information, yielded sig- 
nificantly more hardness ratios above the chosen significance 
threshold. 

4. Images 

For Figs. 2, 3, 4, we smoothed the images and un-vignetted 
exposure maps per band and observation (created as described 
in Sect. 2) with a Gaussian (FWHM of 20") and masked them 
for acceptable detector area. For PN, we subtracted OOT im- 
ages, that were masked and smoothed in the same way as 
the images. We normalized the images with the un-vignetted 
exposure maps to avoid an over-correction of the partly un- 
vignetted detector background. For the EPIC combined images, 
we added the images of the individual cameras scaled accord- 
ing to the background in the individual energy bands. For the 
colour image and the broad band image we added images for 
the individual bands as needed and masked them to a total un- 
vignetted exposure of 10 ks. To better visualize faint structures 
we add contours. We indicate the detected sources by boxes 
and the optical extent of M 33 by the incUnation-corrected D25 
ellipse. 

The XID band (0.2-4.5 keV) image (Fig. 2) gives an 
overview on the sources detected in the XMM-Newton anal- 
ysis. Two of the sources (1 and 7) are bright sources detected 
in areas with low exposure that were masked out for the im- 
ages. In Fig. 3 we rotated the images of the individual band in a 
way that the major axis of the D25 elUpse is oriented vertically. 



We over-plotted the positions of the sources of our catalogue. 
By comparing the different energy bands it is clear that many 
sources only show up in some of the images, a fact that indi- 
cates spectral diversity and is further quantified in the different 
hardness ratios of the sources. This fact can be visualized even 
more clearly in the combined EPIC colour image where we 
coded the 0.2-1.0 keV band in red, 1.0-2.0 keV band in green 
and above 2.0 keV in blue (see Fig. 4). The image is a demon- 
stration of the colourful X-ray sky. SSS, thermal SNRs and 
foreground stars appear red or yellow, XRBs, Crab-like SNRs 
and AGN green to blue. Diffiise emission fills the inner area of 
the disk and the southern spiral arm (red, see also Fig. 3, band 
1-3). A detailed analysis of this emission is outside the scope 
of this paper. 



The images reveal that there is additional emission that has 
not been resolved into individual sources. This may be due to 
really diffuse emission or due to sources that did not reach 
the detection threshold. Examples are, structures in the area 
around the three bright sources X-9a,b,c (277, 297, 3 10) in the 
NE, source X-^ (124) in the NW and source X-6 (155) in the 
SSW. There seems to be a relatively bright source which was 
missed by the detection procedure between sources 235 and 
258. Faint structure appears between sources 55 and 65 and 
south of source 91. Only deeper and/or higher resolution ob- 
servations will allow us to clarify the nature of this emission. 
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Table 4. X-ray source catalogue of the XMM-Newton EPIC M 33 raster observation. 



Src 


Val 


EPIC PN parameters 
Exposure ML^ji 


Count rate* 


Flux* 


HRl 


HR2 


HR3 


HR4 


(1) 


(2) 


(ks) 
(3) 


(4) 


(cts-i) 
(5) 


(erg cm"^ s~') 
(6) 


(7) 


(8) 


(9) 


(10) 


1 


TFF 


2.7 


3.41e+01 


1.18e-02±2.2e-03 


1.55e-14±4.0e-15 


0.44±0.21 


-0.39±0.21 


-0.84±0.34 


*^ * 


2 


TTF 


2.7 


1.32e+02 


3.30e-02±3.7e-03 


4.94e-14±7.2e-15 


0.38±0.12 


-0.14±0.13 


-0.71±0.16 


*^ * 


3 


TTT 


3.8 


1.04e+01 


4.40e-03±1.4e-03 


5.91e-15±2.4e-15 


0.65±0.33 


-0.59±0.32 




*^ * 


4 


TTT 


3.8 


1.46e+01 


4.42e-03±1.4e-03 


1.05e-14±4.1e-15 




0.27±0.34 


-0.32±0.36 




5 


TTT 


8.4 


1.57e+01 






43+n 32 


-0 45+0 31 




-1.00±1.02 


408 


TTT 


4.1 


5.06e+01 




2.14e-14±4.7e-15 




-fl n7+0 14 


-0.54±0.22 


*^ * 






EPIC MOSl parameters 
Exposure ML^ji 


Count rate* 


Flux* 


HRl 


HR2 


HR3 


HR4 






(ks) 
(11) 


(12) 


(cts-i) 
(13) 


(erg cm"^ s"') 
(14) 


(15) 


(16) 


(17) 


(18) 








* 






*^ * 


•± * 










2.8 


1.42e+01 


5.61e-03±1.7e-03 


2.78e-14±9.6e-15 


0.76±0.30 


-0.28±0.3() 


--fi 








9.5 


2.64e-ai 


4.38e-04±4.5e-04 


2.60e-15±3.2e-15 






*± * 








4.9 


1.38e-01 


3.15e-04±4.2e-04 


2.57e-15±4.1e-15 














10.7 


1.65e+00 


5 SOp-04+3 2e-n4 


S 14P-1S+3 2p-15 


*^ * 


+^ * 


*^ * 


*^ * 






4.1 


3.84e+01 


7.34e-03±1.5e-03 


5.23e-14±1.3e-14 


1.00±0.34 


0.23±0.22 


-0.23±0.23 


-0.02±0.39 






EPIC M0S2 parameters 


















Exposure 


MLexi 


Count rate* 


Flux* 


HRl 


HR2 


HR3 


HR4 






(ks) 
(19) 


(20) 


(cts-i) 
(21) 


(erg cm"^ s~') 
(22) 


(23) 


(24) 


(25) 


(26) 






* 


* 




*^ * 


*^ * 












* 


* 










* 


*^ * 






9.4 


1.84e+01 


2.32e-03±6.Ie-04 


1.26e-14±3.8e-15 


0.01±0.26 


-0.72±0.31 




* 






4.7 


4.22e+00 


1.64e-03±7.6e-04 


1.54e-14±7.4e-15 








*^ * 






10.4 


1.83e+00 


9.21e-04±4.5e-04 


5.84e-15±3.1e-15 


*^ * 






*^ * 






4.0 


2.35e+01 


6.09e-03±1.5e-03 


4.43e-14±1.3e-14 


*^ * 


0.17±0.32 


-0.10±0.33 


-0.38±0.58 



Notes: 



*: in XID band (0.2-4.5 keV) 



Fig. 2. Logarithmically-scaled XMM-Newton EPIC low background image of the M 33 medium and thin filter observations combining PN and 
MOSl and MOS2 cameras in the (0.2^.5) keV band. The data are smoothed with a Gaussian of FWHM 20" which corresponds to the average 
point spread function of the combined observations with different off-axis angles. The image is corrected for un- vignetted exposure and masked 
for exposure above 14 ks. Contours are at (0.8, 1.2, 2.0, 3.6) x 10"' ct s"'. Sources from the catalogue, RA, DEC (J2000.0) coordinates, and the 
optical D25 ellipse of M 33 are marked. 



Fig. 3. XMM-Newton EPIC M 33 images: band 1-4 (upper left), band 1 (upper right), band 2 (lower left), and band 3 (lower right). 



5. Cross-correlation with other M 33 X-ray 
catalogues 

In this section we discuss the cross-correlation of the XMM- 
Newton detected sources with sources reported in earher X- 
ray catalogues. All correlations (together with other X-ray in- 
formation like variability, reported in these catalogues, or ex- 
tent, detected in this work) are indicated in the XID column of 
Table 3 (col. 13). 

From the 17 Einstein sources reported in Trinchieri et al. 
(1988) we detected all but the relatively bright source X-12 in 
the northern disk and the less than 3cr detection X-15 which 
was located about 90" to the East of the nucleus. Both of 
these sources were also not detected by ROSAT (Schulman 
& Bregman 1995; Long et al. 1996, HPOl). While one may 
question the original existence of X-15, X-12 was detected by 
Einstein with high significance and has to be strongly time vari- 
able. 



In the correlation with ROSAT detected sources, we con- 
centrated on the catalogue derived by HPOl from all archival 
observations of the field, which contains 184 sources within 
50' of the nucleus. It re-detects all 37 sources detected in the 
ROSAT PSPC analysis of Long et al. (1996) and 21 of the 27 
sources from the first ROSAT HRl analysis of the central area 
of M 33 by Schulman & Bregman (1995). We searched for 
these low signal-to-noise sources (numbers 18, 19, 20, 23, 24, 
and 27) without success in the more sensitive XMM-Newton 
data. The sources therefore most Ukely are not just time vari- 
able, but rather false detections. 

In Table 5, we classified sources of HPOl who analyzed 
all archival ROSAT PSPC and HRl observations of the M 33 
field: 19 sources are positioned outside the field covered by the 
XMM-Newton observations, three are detected within the ex- 
tended emission of XMM-Newton source 125, 5 are resolved 
into multiple components in the XMM-Newton catalogue. We 
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Fig. 3. (continued) XMM-Newton EPIC M 33 images: band 4 (left) and band 5 (right). The images are smoothed with a Gaussian of FWHM 
20", corrected for un-vignetted exposure and masked for exposure above 14 ks. Contour levels are (8, 12, 20, 36) x 10"* ct s"* for the merged 
band 1-4, (1.5,2,3,5) x IQ-" ct s"' for the band 1, (1.7,2.2,3.2,5.2) x 10"* ct s"' for the band 2, (3,4,6, 10) x IQ-*^ ct s"' for the band 3, 
(2.5, 3.5, 5.5, 9.5) x 10 * ct s ' for the band 4, and (5, 6, 8, 12) x 10 ct s ' for the band 5 images. The inclination-corrected optical D25 ellipse 
of M 33 , the sky orientation, and the sources from the catalogue - including source numbers in the (4.5-12) keV band image - are indicated. 



Fig. 4. Logarithmically-scaled, three-colour 
XMM-Newton EPIC low background im- 
age of the M 33 medium and thin fil- 
ter observations combining PN and MOS 1 
and M0S2 cameras. Red, green and blue 
show respectively the (0.2-1.0) keV, (1.0- 
2.0) keV and (2.0-12.0) keV bands. The 
data in each energy band have been 
smoothed with a Gaussian of FWHM 20" 
which corresponds to the average point 
spread function of the combined observa- 
tions with different off-axis angles. The im- 
age has been corrected for un-vignetted ex- 
posure and masked for exposure above 14 
ks. The image scale and the inclination- 
corrected optical D25 ellipse of M 33 are 
marked. 



Table 5. Summary of cross-correlation with the M 33 ROSAT cata- 
logue of HPOl. ROSAT sources, not detected in our XMM-Newton 
analysis, are arranged according to their ROSAT detection likeUhood 
LH. 



classification 



HPOl source numbers 



outside XMM-Newton field 

within extended source 125 

resolved by XMM-Newton 

not detected: LH<9 
(28 of 34 sources) 

not detected: 9<LH<10 
(13 of 19 sources) 

not detected: LH>10 



1,2,3,4,5,6,7,10,11,46,167,171, 
176,177,180,181,182,183,184 

48,50,51 

22,29,59,121,172 

14,16,24,35,37,38,65,73,79,85,88, 

89,94,105,107,108,116,127,131, 

135,138,149,153,155,159,166,175 

23,36,43,47,56,72,90,99,100,101, 
113,125,142 

9(LH=29),30(59),41(44),53(13), 
68(42),69(1 1),82(27),83(14), 
93(198,near nucleus),97(ll), 
98(21,near nucleus), 119(14,near 
nucleus),128(18),133(20) 



did not detect 41 out of 51 sources with a ROSAT detec- 
tion likelihood below 10. All these sources originate from the 
ROSAT HRI detections list which, following the experience at 
the time, was generated with a likelihood limit of 8 for detec- 
tion compared to the limit of 10 used for PSPC sources. If still 
at the ROSAT brightness, these HRI sources should have been 
detected by the deeper XMM-Newton observations. While, in 
principle, they could all have dimmed by such an amount that 
we did not detect them, this seems rather unlikely and prob- 
ably most of these sources were spurious detections. The sit- 
uation is certainly different for the HPOl sources with likeli- 



hood above 10. Many of these sources were detected with both 
ROSAT detectors and the non-detection with XMM-Newton 
indicates time variability. HPOl source 30 was already found 
to vary by more than a factor of 7.5 during the ROSAT obser- 
vations, two other sources (69 and 82) were classified by HPOl 
as SSS candidates, which as a class are known to show strong 
time variability, and another two were classified as foreground 
stars (9 and 133), which might have been flaring during the 
ROSAT observations. The HPOl sources close to the nucleus 
of M33 (93, 98, and 119) may have been missed in our XMM- 
Newton analysis due to contamination from the bright nuclear 
source X-8 or could be spurious detections due to structures 
in the ROSAT PSF. One other HPOl source (128) is clearly 
visible in the XMM-Newton band 3 and 4 images (see Fig. 3) 
between XMM-Newton sources 235 and 258, however it was 
not picked up by the detection algorithms. Besides these spe- 
cial cases, there are 92 clear and four questionable (214, 256, 
258, 338) XMM-Newton and HPOl correlations listed in col. 
13 of Table 3. 



6. Classes of point-like X-ray sources detected in 
the direction of M 33 

To identify the M 33 X-ray sources we searched for corre- 
lations around the X-ray source positions (within a radius of 
3 X {(Tstat + CTsyst) in the SIMBAD and NED archives and within 
several catalogues. In cols. 14 to 19 of Table 3, we give extrac- 
tion information from the USNO-Bl catalogue (name, number 
of objects within search area, distance, B2, R2 and 12 magni- 
tude of the brightest object). We omitted information, when the 
digitization was fooled by extended emission within Hn re- 
gions (mostly only some colours available and rather bright). 
To improve on the reliability of identifications we used the B 
and R magnitudes to calculate the ratio of X-ray- to-optical flux. 
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Table 6. Summary of identifications and classifications. 





Source type 


Sclcclion criteria 


idcnlilicd classified 



fg Star log(;^^) < -1.0 and HR2 < 0.3 and HR3 < -0.4 or not defined 5 30 

AGN Radio source and not classification as SNR from HR2 12 

GAL optical id with galaxy and HR2 < 0.0 11 

SSS HRl < -0.2, HR2 - EHR2 < -0.99 or HR2 not defined, HR3, HR4 not defined 5 

SNR HRl > 0. 1 and HR2 < -0.4 and not a fg Star 21 23 

XRB optical identification or X-ray variability 2 

hard HR2-EHR2 > -0.2 or only HR3 and HR4 defined, and no other classification 267 



given by log(i^) = log(fx) + (mB2 + mR2)/(2 x 2.5) + 5.37, fol- 
lowing Maccacaro et al. (1988, see col. 20). 

The catalogued X-ray sources are identified or classified 
based on properties in X-rays (hardness ratios (HR), variability, 
extent) and of correlated objects in other wavelength regimes 
(Table 3, cols. 21, 22). The criteria used are summarized in 
Table 6. We counted a source as identified, if at least two cri- 
teria secure the identification. Otherwise, we only counted a 
source as classified (indicated by pointed brackets). 

We plotted X-ray colour/colour diagrams based on the HRs 
(see Fig. 5). Unclassified sources are only plotted if the error in 
both contributing HRs is below 0.2. To identify areas of specific 
source classes in the plots, we over-plotted colours of sources, 
derived from measured XMM-Newton spectra and model sim- 
ulations (see discussion in subsections below). 

Identification and classification criteria and results are dis- 
cussed in detail in the subsections on individual source classes 
below. Many foreground stars, SSS and SNRs can be classi- 
fied or identified. However, besides a few clearly identified 
XRBs and AGN, and SNR candidates with known positions 
at other wavelengths, we have no clear hardness ratio criteria 
(see below) to select XRBs, Crab-like SNR or AGN. They are 
all "hard" sources and we therefore introduced a class <hard> 
for sources with HR2 minus HR2-error greater than -0.2 (see 
Tab. 6). 

43 sources remained unidentified or without classification. 
6.1. Foreground Stars 

Foreground stars are a class of X-ray sources that is homo- 
geneously distributed over the FOV of the M 33 survey. The 
good positioning of XMM-Newton and the available catalogue 
USNO-Bl allowed us to effectively select this type of sources. 
We identified five sources with stars of known type in the 
SIMBAD data base. Their X-ray fluxes were in the range ex- 
pected for these sources. In addition, we classified 30 sources 
as foreground stars (log(^) < -1, see Maccacaro et al. (1988) 
and in addition HR2 < 0.3 and HR3 < -0.4). For several of the 
star candidates we estimated the type from the colours in the 
USNO-Bl catalogue using the steUar spectral flux library from 
Pickles (1998). For ten candidates we could not determine the 
type from the colours (28, 54, 77, 122, 128, 262, 308, 337, 371, 
393). This may indicate that these sources are not isolated stars 
but more complicated systems or even, in some cases, galaxies. 



This expectation could, for instance, be confirmed for source 
262, which was resolved into two sources in DSS2 images. 

A special case was source 196, only ~ 2' WSW 
of the nucleus. The source correlated with 2MASS 
J01334186+3038491, and was catalogued by Ivanov et al. 
(1993) as red star 184. The X-ray HRs, optical colour and 
log(7^) point towards a foreground star of type M6 as the 
identification of this source. 

From the HRs, sources 7 and 10 should be foreground stars. 
However, there was no candidate star within the error box. As 
these sources were at the border of the field, the systematic er- 
ror may be underestimated and nearby stars may be candidates 
for the optical identification. 

For source 240, the HRs indicated a foreground star, how- 
ever log(i&-) - -0.7, larger than the assumed limit for clas- 

I op I 

sification. The optical colours would indicate a MIO III star. 
Also the other candidate stars of type late M have a rather high 
log(^). 

^opt 

For several sources, log(5^^) pointed towards a stellar iden- 

Topi 

tification (8, 38, 72, 176, 206, 209, 226, 234). However, only 
for three sources the HRs including errors were within the lim- 
its we assume for foreground stars (72, 176, 234). 

6.2. Galaxies and AGN 

Already after the Einstein observations, the X-ray bright source 
X-9c (297) was identified by Christian & Schommer (1982) 
with an elliptical galaxy at z=0.03. The XMM-Newton X-ray 
luminosity of 1.35x10"*^ erg s^^ (0.2^.5 keV) was typical and 
the HRs were compatible with an elliptical galaxy spectrum 
(see HR plots). 

From SIMBAD and NED searches, we found an additional 
galaxy correlation (nameley source 358 - LEDA5899). 

There were no correlations with AGN with known red- 
shift within SIMBAD or NED. We therefore only classifed 
12 sources as AGN based on SIMBAD, NED, and other ra- 
dio source correlations (NVSS, Gordon et al. (1999)) with the 
additional condition of being a <hard> X-ray source. A final 
decision will only be possible using optical spectra. In Fig. 5, 
we included typical HRs for an AGN spectrum (power law with 
photon index of 1 .7 assuming galactic foreground absorption 
(6x10'^'^ cm"^, Al) and absorption through the densest areas of 
M 33 (Newton 1980, 4.8x10^1 cm'^, A2). 
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Fig. 5. Hardness ratios of sources detected by XMM-Newton EPIC. 

Shown as dots are only values with errors smaller than 0.20 on both 
HR(i) and HR(i+l). Foreground stars and candidate are marked as big 
and small stars, galaxies and candidates as big and small x, AGN can- 
didates as crosses, SSS candidates as triangles, SNR and candidates as 
big and small hexagons, XRBs as squares. In addition, we mark posi- 
tions derived from measured XMM-Newton EPIC spectra and models 
for SSS (SI to S4) as filled triangles, low mass XRBs (LI and L2) as 
filled squares, SNRs (N132D as Nl, IE 0102.2-7219 as N2, N157B 
as N3, Crab spectra as CI and C2) as filled hexagons, AGN (Al and 
A2) as asterix. 



Sources 150 and 373 could be classified as SNRs, as well 
as AGN. They were identified as radio sources within the 
inclination-corrected M 33 D25 ellipse, without optical iden- 
tification as SNR by Gordon et al. (1999) and the radio spectral 
index did not clearly classify them as AGN. 

6.3. Super-Soft Sources (SSS) 

SSS show black body spectra with temperatures below 50 eV, 
they radiate close to the Eddington luminosity of a 1 Mo object 
and they are believed to be white dwarf systems steadily burn- 
ing hydrogen at the surface. They were identified as a class of 
X-ray sources by ROSAT and are often observed as transient X- 
ray sources (see Greiner 2000, and references therein). In the 
catalogue, SSS were only classified using their HRs. To guide 
the classification, we calculated SSS HRs assuming a 25 and 50 
keV black body model spectrum assuming galactic foreground 
absorption (SI and S3) and absorption through the densest ar- 
eas of M 33 (S2 and S4). This led to the selection criteria in 
Table 6 and the classification of five M 33 SSS (52, 109, 223, 
247, 332). They were detected with absorbed luminosities of 
(l-8)xl035 5-1 ti^e 0.2-0.5 keV band. Only three of the 
sources (and the model HRs) could be plotted in Fig. 5. For the 
other two, HR2 was not determined. 

6.4. Supernova remnants (SNR) 

SNRs can be identified as sources where thermal components 
dominate the X-ray spectrum below 2 keV (examples are 
N132D in the Large Magellanic Cloud and IE 0102.2-7219 
in the Small Magellanic Cloud) and as so called "plerions" 
with power law spectra (examples are the Crab nebula and 
N157B in the Large Magellanic Cloud). To guide the classifi- 
cation we calculated HRs from archival XMM-Newton spectra 
of these SNRs. Spectra of N132D (Nl), IE 0102.2-7219 (N2), 
and N157B (N3) could be directly compared to M 33 SNRs, as 
they are seen through comparable foreground absorption. For 
the Crab nebula spectrum, we assumed galactic foreground ab- 
sorption (CI) and absorption through the densest areas of M 33 
(C2). It is clear from Fig. 5, that "thermal" SNRs are located 
in areas of the X-ray colour/colour plots that only overlap with 
foreground stars. If we assumed that we have identified all fore- 
ground star candidates from the optical correlation and inspec- 
tion of the optical images, the remaining sources could be clas- 
sified as SNRs with the criteria given in Table 6. 

Extensive searches for SNRs in M 33 were performed in the 
optical (see Gordon et al. 1998, and references therein) and ra- 
dio (see Gordon et al. 1999, and references therein). Many X- 
ray sources in the catalogue correlated with optical and/or radio 
SNRs (see Tab. 3). We counted such sources as SNR identifi- 
cations (21), if in addition they fulfilled the HR criteria (see 
above). We counted sources as classified SNRs (23), if they ei- 
ther just fulfilled the HR criteria (16) and were no identified 
or classified as foreground star, or if they correlated with opti- 
cal/radio SNRs and did not fulfill the HR criteria (7). 

From these classified SNRs, sources 20 fulfilled the HR cri- 
teria, but was positioned outside the incUnation-corrected op- 
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Fig. 6. XMM-Newton EPIC 0.2-2 keV band contours of 
(6,8,14,32,86)xlO"^ ct s"' overlayed on an optical image of 
M 33 extracted from the digitized sky survey (DSS2) red plates. The 
image scale is marked. 

tical D25 ellipse, while sources 101 and 266 were located be- 
tween two optical SNR candidates and may represent the com- 
bined emission of both. As already mentioned in the subsection 
on AGN, we classify sources 150 and 373 as SNR as well as 
AGN due to the radio correlation. 

A special case was source 227, which correlated with ra- 
dio source 126 of Gordon et al. (1999), who classified it as a 
background source with no optical identification. The XMM- 
Newton HRs clearly indicated a SNR nature and we therefore 
classified the source as a SNR. 

The giant H 11 regions BCLMP 290 and NGC 604 were de- 
tected as extended X-ray sources (125 and 299). While the HRs 
of source 299 indicated a spectrum similar to thermal SNRs 
(several SNRs would be necessary to explain the X-ray emis- 
sion), the HRs of source 125 indicated a hard spectrum which 
may originate from several point sources (XRBs) or Crab-like 
SNRs. In the catalogue we only classified source 299 as a SNR. 

Four additional sources in the catalogue (137, 154, 162, 
280) had HRs compatible with the selection criteria in Table 6, 
within the errors. They are good candidates for further SNR 
investigations in M 33. 

6.5. X-ray binaries (XRB) 

As already mentioned in the introduction to this section, ex- 
pected spectra of XRBs are similar to AGN and Crab-like 
SNRs. To guide the classification, we calculated HRs as ex- 
pected for low mass XRBs (5 keV thermal Bremsstrahlung 
spectrum assuming galactic foreground absorption (LI) and ab- 
sorption through the densest areas of M 33 (L2)), which show 
soft spectra for XRBs. As can be seen in Fig. 5 these different 
source classes did not separate. 

Therefore we only listed the black hole XRB X-8 close to 
the nucleus (211) and the 3.4 d eclipsing HMXB X-7 (171) as 
identified X-ray binaries. 

SIMBAD searches indicated 16 more good XRB candi- 
dates from positional coincidences with optical stars of differ- 
ent type in M 33 (see Table 3). Source 110 correlated with an 
emission fine object and is a good candidate for a Be XRB, 
source 174 correlated with an echpsing binary candidate (398 
d period), and sources 194, 197, and 323 with variable stars. 
Of historical interest is the optical candidate for source 323, 
VHK6. This source was discussed as one of four irregular vari- 
ables in the paper by Hubble (1926), in which he proposed the 
spiral nebula M 33 as a stellar system outside the Milky Way. 

AH the candidates mentioned above need further X-ray 
work to confirm their XRB nature. 

7. Discussion 

Besides the large number of point-like sources, there was a 
strong component of unresolved X-ray emission in the energy 



band below 2 keV. The overlay of X-ray on the DSS2 red im- 
age (Fig. 6) clearly showed the strong correlation with the inner 
disk and the southern and NW spiral arm. While some of this 
emission certainly is due to unresolved point sources, there is a 
significant contribution from hot gas in the disk which is more 
concentrated in the spiral arms and H 11 regions like NGC 604 
12' NE and BCLMP 290 15' NW of the nucleus. A detailed 
analysis of this emission will be presented in a separate paper 

In the XMM-Newton source catalogue we presented 408 
sources in the M 33 field, 35 of which could be identified or 
classified as foreground stars and 14 as background AGN or 
galaxies. As sources in M 33 we identified or classified 5 SSS, 
44 SNRs, and 2 XRBs. 267 sources were classified as hard and 
may be either XRBs or Crab-Uke SNRs in M 33 or background 
AGN. There remained 43 sources unclassified. 

We can estimate the contribution of background sources to 
the source catalogue extrapolating from X-ray deep field num- 
ber counts in the 0.5-2.0 keV and 2.0-8.0 keV bands (see e.g. 
Brandt et al. 2001; Rosati et al. 2002). We used rough sensi- 
tivity maps for our survey created by the SAS task esensmap 
for the 0.5-2.0 keV and 2.0^.5 keV bands and an estimate for 
additional absorption in the hne-of-sight. One expects about 
175 to 210 background sources detected in our field-of-view, 
if we assume an overlap of sources detected in the hard and 
soft band of 75% or only 50%. Fourteen of the background 
sources have already been identified (see Sect. 6.2). Taking 
these numbers into account, up to 50% of the hard or unclassi- 
fied sources should be sources within M 33 (mostly XRBs or 
Crab-like SNRs). 

We compared the XMM-Newton source catalogue with the 
ROSAT catalogue of HPOl. While most of the bright ROSAT 
sources were also detected by XMM-Newton many of the 
sources, only detected with the ROSAT HRI with a detection 
likelihood below 10, were not detected by XMM-Newton and 
seem to be spurious detections. From the seven ROSAT sources 
classified as SSS, we did not detect four As expected for this 
source class, these might be transients, that were not active dur- 
ing the XMM-Newton observations (HPOl sources 43, 69, 76, 
82). The remaining three seem to be false identifications: HPOl 
source 61 is now classified by XMM-Newton as a SNR, 84 as a 
foreground star and 124 as a hard source. Two ROSAT correla- 
tions with radio sources are no longer valid with the improved 
XMM-Newton positions (HPOl sources 63 and 147). 

HPOl identified and classified 17 SNRs. The XMM- 
Newton source catalogue marks 21 identifications and 23 clas- 
sifications. While the HR criteria seem to be rather robust 
for thermal SNRs, the identification of Crab-like remnants is 
doubtful and can only be secured by measurement of source ex- 
tent. This may be possible with XMM-Newton EPIC for very 
extended SNRs (optical extent of SNR in M 33 vary from 9" to 
126" according to Gordon et al. (1998)). Most of the sources, 
however, will need dedicated Chandra on-axis observations for 
a final decision on the extent. The Chandra count rates will 
not be sufficient to get good spectra and a combined XMM- 
Newton and Chandra analysis will be important for this class 
of sources. 

The importance of combined XMM-Newton and Chandra 
data for M 33 XRB research was demonstrated in Pietsch et al. 
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(2004). This type of analysis extended to all brighter sources 
will be the subject of a future paper. 
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